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Abstract 

Correlated Basis Function theory and Fermi Hypernetted Chain technique are 
extended to study medium-heavy, doubly closed shell nuclei in j-j coupling scheme, 
with different single particle wave functions for protons and neutrons and isospin de- 
pendent two-body correlations. Central semirealistic interactions are used. Ground 
state energies, one-body densities, distribution functions and momentum distribu- 
tions are calculated for ^^C, ^^O, ^°Ca, ^^Ca and ^osp)^ nuclei. The values of the 
ground state energies provided by isospin dependent correlations are lower than 
those obtained with isospin independent correlations. In finite nuclear systems, the 
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two-body Euler equations provide correlation functions variationally more effective 
than those obtained with the same technique in infinite nuclear matter. 



1 Introduction 

The description of the properties of all nuclear systems, from deuteron to nuclear matter, 
with a hamiltonian containing realistic two- and three-body potentials is one of the 
goals of non relativistic nuclear many-body theory [|l|. We call realistic those two-body 
potentials which reproduce the ground state properties of the deuteron and the nucleon- 
nucleon scattering data. Three-body forces are generally required to provide a good 
description of three-body nuclei. 

The techniques to solve the Schrodinger equation in few-body systems with these 
potentials have reached a high degree of sophistication and have produced excellent re- 
sults [0]. Accurate and reliable technologies have also been developed for infinite nuclear 
systems, and, nowadays, their results are quite satisfactory P, 

The situation for medium and heavy nuclei is still troublesome. The extension to these 
nuclei of the Monte Carlo techniques used in few-body systems is hindered by computa- 
tional difficulties related to the relatively high number of particles involved. In a nucleus 
like ^'^Ca the number of spin isospin configurations to be sampled is of the order of the 
Avogadro's number. On the opposite side, the techniques which made possible accurate 
numerical studies in infinite nuclear and neutron matters cannot be straightfordwally 
applied to finite nuclear systems because they do not enjoy translational invariance. 

In two previous papers p, |^ we succeeded in extending Correlated Basis Function 
(CBF) theory and Fermi Hypernetted Chain (FHNC) cluster summation technique 
(successfully used to describe nuclear matter properties) to deal with finite Fermi systems. 
In these works the FHNC theory has been applied to study model, doubly magic nuclei 
with the same number of protons and neutrons, switching off the Coulomb interaction 
and using single particle bases in the 1-s coupling scheme. 

This scheme can be adequate to describe closed shell nuclei up to ^°Ca, but it breaks 
down for heavier nuclei, where the number of neutrons is larger than that of protons. 
Moreover, while the shell closure for A < 40 corresponds to that of 1-s coupling scheme, 
for heavier nuclei it is necessary to use a basis with spin-orbit splitting to reproduce the 
correct sequence of magic numbers. 

In the present paper we make a step forward towards a fully microscopic description 
of all doubly closed shell nuclei. We extend the FHNC formalism to differentiate protons 
and neutrons and, in addition, to consider single particle bases in a j-j coupling. 

In the next section we shall present the modifications of the FHNC equations needed 
to describe the above nuclear systems. 

Nuclei are described as a mixture of neutrons and protons, with possibly different 
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populations. This is reflected by the presence of distribution functions depending on the 
isospin of the reference nucleons. In a single particle basis generated by a spin-orbit term, 
the third components of the spin is no longer a good quantum number to describe the 
single-particle wave functions. For this reason, a new type of statistical link appears in 
the FHNC equations. 

We shall discuss the influence of the aforementioned FHNC modifications on the cal- 
culation of the nuclei ground state energy in section 3. 

The results of the calculations of binding energies, density and momentum distribu- 
tions, for various doubly closed shell nuclei are presented in section 4. These calculations 
have been performed using central nucleon-nucleon interactions (^4) with spin and isospin 
dependence and without tensor components. The correlation functions have isospin de- 
pendence. Since these interactions and correlations are not yet fully realistic, the obtained 
results cannot be immediately compared with the experimental data. These are, however, 
the first, true microscopic calculations in such heavy nucleus as ^ospi^ they represent 
a necessary step towards an ab initio description of heavy nuclear systems starting from 
a realistic nuclear hamiltonian. 



2 Extension of the FHNC method 

In the framework of CBF theory and along the lines of the formalism presented in |^ , the 
ground state of a nucleus with Z protons and neutrons {A = Z + N) is described by 
the correlated wave-function: 

vl>(l,...,/l) = F(l,...,/l)$^,z(l,...,A) (1) 

where F is an A-body correlation operator we shall specify later. The function $ n^z is a 
Slater determinant of single particle (s.p.) wave functions (a = p,n) generated by 

the s.p. hamiltonian 

^"p(0 = -^V? + f/"(0- (2) 

A proper description of the sequence of magic numbers is achieved only if the s.p. hamil- 
tonian contains a spin-orbit interaction. In this case, the s.p. wave functions are classified 
in terms of the total angular momentum j and are eigenvectors of and jz- We express 
them as: ^ 

where Xs are the spin wavefuntions, Yj^ the spherical harmonics and (//il/2s | jm) the 
Clebsch-Gordan coefficients. 

The FHNC equations can be written in terms of the one-body densities and the two- 
body distribution functions : 

P?(r) = (^*E^(r-r.)n"*) (4) 

k=l 
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pg(r,r') = J: 6iv-r,)P^6ir'-ri)Pl^Ol,^) g=l,..4, (5) 

where is the projection operator on the a = p,n state of the /c-nucleon and we have 
defined 

dr meaning integration over the spatial coordinates as well as sum over the spins. The 
index q = 1, .A labels the operatorial component of p2.g- This dependence is the same as 
that of the V4 potentials adopted, 



v,il,2)= Y: v^'KnM, (7) 

9=1,4 



with = 1, 0^2 = CTi ■ <^2, = "^1 ■ "^2 and Of-^ = cri ■ (T2T1 ■ T2. 
The uncorrelated density matrices are given by: 

p^(ri,r2) = Y.<Priri)<Pt{r2) 

k 

= Epf"(ri,r2)xl(l)x.'(2), 

s,s' 

and their spin dependent parts are: 

n,l,j 



J2 (^^^^ I Jm){lf^'^s' I jm)Yi*^{fi)Yi^'if2). (9) 

fj,,fj.',m 

While in the \Ps'^lzSz > representation Pq^'" is diagonal in the spin variables in the 
|^2_g2j2j^ > one it depends on the third components of spin (s, s'). The same dependence 
is also present in the correlated density matrices, 

vi/*(l,2,...,A)vl/(l',2,...,A)xI,(l')xI(l'), (10) 

where a sum over all the spin coordinates is implied. 

We found convenient to consider separately the uncorrelated density (|) for pairs of 
particles having either the same or opposite third spin component. For this reason we 
define 



Plp{r^,r,) = -5:(2j + lX,^.(rO<,,(r2)PKcos^i2), (11) 



1 
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where Pi{x) is the Legendre polynomial of 1-th degree and Pl{x) is its first derivative 
with respect to x. The function defined in Eg . ([TT|) corresponds to pairs of nucleons with 

parallel spins, s = s' ( ppp = pg ^ = Po ^ ^ ) ■ This is the only statistical correlation 

appearing in 1-s coupling. The function defined in Eg. (]T2|) , = Po = — Po is 
instead a new statistical link between nucleons with antiparallel spins, s = —s', due to 
the j-j coupling. This new function is antisymmetric under the exchange of s with s' and 
of the spatial coordinates. In all the observables whose mean value we have calculated 
(energy, density, momentum distribution), we found that the contribution of pg p is much 
larger than that of Pq^. As a matter of fact, if the nucleus is closed in 1-s coupling (as for 
^^O and ^°Ca) and no spin-orbit term is included in the s.p. potential, the contribution 
of Pq,a vanishes. 

The uncorrelated statistical functions appear in the cluster expansion of the two- 
body distribution functions and of the one-body density matrices forming closed exchange 
loops. In 1-s coupling, all the nucleons involved in the loops must have the same third 
spin component. In j-j coupling this is no longer true, as there may be statistical links 
between particles with opposite third components. In general, the contribution from a 
N-particle exchange loop is given by: 

(-1)^-V?(ri, r2)p?(r2, v,) . . . p^(r^, n) = 

i-lf-' pr"(ri,r2)pr"(r2,r3)...pr^"(r^,ri). 

S1,...,SM 

3 Energy, density and momentum distribution ex- 
pectation values 

In the previous section we have discussed the modifications induced by j-j coupling and 
Z ^ N on the uncorrelated density matrix. In this section we will shortly discuss how the 
FHNC eguations change and how the expectation values of energy, density and momentum 
distribution are calculated. 

The structure of the FHNC eguations depends on the adopted correlation function. We 
use an isospin dependent correlation, since we want to distinguish protons from neutrons, 
and we choose the following form for F, which allows us to consider different correlations 
for different pairs of nucleons: 

F(l,...,A)= n Uin,)PrP,^. (13) 

l=i<j 

The form of the correlation operator (|13|), and specifically the fact that the involved 
operators commute among themselves, allows us to use, in the cluster expansion, the 
diagrammatic rules presented in Ref. [^, with only slight modifications due to the j-j 
coupling scheme. Since the system is a mixture of different fermions, the correlated density 
matrices for pp, pn and nn are different. Furthermore, we have to consider the presence 
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in the cc-FHNC chains of the new statistical hnk Pg^. The new FHNC expressions of 
the two-body distribution functions and of the one-body density matrices are given in 
the Appendices A and B, respectively. 

The calculation of the expectation value of the energy has been performed along the 
same lines followed in Ref. 0. 

The evaluation of the linetic energy has been done using the Jackson-Feenberg. This 
allows us for eliminating terms of the form (Vi-F)(Vi$o); involving three-body operators 
expectation values P). Following the same notation of Ref.0, we express the kinetic 
energy as: 

< T >= TjF = Tp + T^'^ + rf + . (14) 

In the previous equation we have indicated with Tp those terms where the kinetic 
energy operator acts on the correlation operator F only. Their contribution is given by: 



h_ 

n 

a, 13 

where t[f] has been defined as: 



TF = -^Y.j d'nd'r,pf,ir^, r,)t[Mn,)], (15) 



t[f{r)] ^ y^{/(OrW + Imfir) - r(r)}. (16) 

The remaining j'j""^'^'^) components are the contributions of those terms containing the 

operator acting on the uncorrelated s.p. wave functions. T^^'' is the sum of the cluster 
diagrams in which the external point 1, argument of the one-body kinetic energy operator 
Ti, is not involved in any exchange (see Eq. 2.24 in [^). We obtain the expression: 

tS,'^ = 5 / ^'^iP?i(ri)C(ri), (17) 

where p^^ is defined as: 

PTi(ri) = EC*(ri)V?0r(ri) " E Vi0r(ri) ■ Vi0r(ri), (18) 

i i 

and is defined in appendix A. 

The other two parts (T^^'^-*) contain diagrams where the point 1 belongs to a 2-body 
(T^^) or to a n{> 2)-body (T^^"*) exchange loop. We find for T^^ the expression: 

= i^E/ d'n£r,C{ri)pUrur2) [gZiri,v,)Cir,) - 1] , (19) 

where 

PT2{ri, rs) = Po,p(ri, r2)Vlplp{ri, rs) - Vip[^_p(ri, rs) ■ Vip[^_p(ri, rs) + 

p^^in, r^MplAri, r2) - Vip^,^(ri, r^) ■ Vip^,^(ri, r^). (20) 
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All the FHNC quantities introduced in the above equations are defined in appendix C. 
The evaluation of T^^-* requires the knowledge of the three-body distribution functions, 

but its leading term (Tj^'^"*) is a function of two-body dressed FHNC quantities 0. T^^'^"* 
is given by: 
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Tf''^ = '^11 j [pT3,p(ri, r2)//rc,p(ri, r^) + p?3,A(ri, r2)ifrc,A(ri, r^)] , (21) 
where 

^cV(ri,r2) = C(ri){C(r2)[(^7r(ri,r2) - l)iV° ^(ri, r^) + iVi$(ri, r^)] 

+ (C(r2) - l)ivS(ri, r^) + Cir2)gZin, r2)E'^^^^{v„ v,) }, (22) 

with X = P,A. The remaining term xj^'^^ contains a three-body operator and it is 
known to be negligible, both in nuclear matter and N=Z nuclei ||^. Therefore, it has been 
neglected. 

The center of mass kinetic energy, T^m, has to be subtracted from < H > to get the 
energy mean value E. T^m is given by: 

Ten. = E / ^'^1 (pTi(ri) - / d'r.pUri, r^)) • (23) 

The expressions of the p^^ functions are given in appendix C. 
The potential energy < V > may be written as: 

<^>=l E / d'^id'r, E v^^\n2)pt,{r„ r,). (24) 

a,l3=p,n 9=1,4 

The two-body distribution functions can be decomposed into a direct part, summing 
those contributions where the particles 1 and 2 do no belong to the same exchange loop, 
and two exchange parts, containing diagrams where the two nucleons are involved in the 
same loop with identical {pexc,p) or opposite {pexc,A) third spin components: 

P2,g(ri, ra) = <^Pdf (^i, r2) + bf p'^^^^pivi, rs) + Pexc,A(ri, rz). (25) 
The a, b and c factors are given in Tab.l. The expressions for p^j^ and Pexc,x are: 

Pf(ri,r2) = e,"(ri)(e,"(r2)^?r(ri,r2)+ef(r2)^?f(ri,r2)) + 
C(ri)d(r2)^7et(ri,r2) + 



C(ri)e(r2)^7r(ri,r2 



iVef(ri,r2) + &,(ri,r2) + 



W/(ri, r2) + i^f (ri, r2))(iV,f (r,, r^) + Ci^(ri, r^)) 



(26) 
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Pefc,x(ri,r2) = C(ri)^f(r2)^dd(ri,r2) 



(<^(ri, r2) + <^(ri, r^) - pj^(ri, r^)) 



(27) 



We consider also the contribution of the Coulomb interaction acting between the 
protons. Its expectation value is: 



2 J ri2 ' 



The mean value of the energy is then given hj E = Tjp+ < V > + < Vc > —Tcm- 
Since the one-body density matrices depend on the third components of the spin 
(Eg. (pup), also the momentum distributions have this dependence : 

n^^'"(A;) = jj c/VirfVre^''-^'"!-"!' V'"(ri, ry). (29) 

The spin averaged momentum distribution, summed over all the possible spin components, 
is: 

n"(fc) = ^n^^'"(A;), (30) 

s,s' 

and, in the following, we shall always refer to the momentum distribution as defined by 
Eg. (pOD . The FHNC derivation of the one-body density matrices are given in appendix 
B. 

We conclude this section by listing some sum rules (SR) that must be satisfied by the 
densities, the two-body distribution functions and the momentum distributions. They are 
particularly relevant since provide information about the accuracy of the approximations 
used to solve the FHNC eguations. 

For the densities and the distribution functions, the following SR's hold: 



Sp 



Sn 



'-'pp 



Snp(pn) 



i / dhpliv) = 1, (31) 

1 J d\p-{v) = 1, (32) 

^^^-^ J rf\irfV2pfi(ri, r2) = 1, (33) 

^-Jd'nd'r,pZ^'-\r,,r,) = l, (34) 

J d'nd'r^pZiru r2) = 1. (35) 



A^(A^ 



In nuclei, where the levels with the same value of the orbital angular momentum are 
saturated (as ^^O and ^°Ca), the spin-SR S^, 
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S^ = ^J d'nd'r [p^^2(ri, rs) + p^slri, r^)], (36) 

has to satisfy the condition S^^ = —1. This is no longer true for the nuclei ^^C, ^^Ca 
and ^°'^Pb. However, as the correlations do not contain spin flip terms, 5*0- in the correlated 
nucleus must have the same value as in the uncorrelated one (5'^°''^ = S^"''^). 

The momentum distributions must obey the zeroth momentum SR: 

MD^ = — i— / £kn'^{k) = Xo,, (37) 



(27r) 



where Xa is the proton {xp = Z/A) or neutron fraction (x„ = N/A), and the second 
momentum SR: 

/ £kk'^[ZnP{k) + Nn^'ik)] =<T> . (38) 



MD 



2m{2n 

4 Results 

We have applied the theoretical framework presented in the previous sections to the study 
of the ground state properties of the ^^C, ^^O, ^°Ca, ^^Ca and ^''^Pb nuclei. 

Our calculations have been done using two different kinds of correlation functions. 
One has a gaussian behaviour: 

/f^(r) = 1 - A„^exp(-5„^r2), (39) 

where Aajj and Bajj are taken as variational parameters, fixed by minimizing the FHNC 
energy. 

The second correlation function is obtained by the minimization of the energy at the 
second order of the cluster expansion, < H2 >, with 



<H2> = E j dh,i^pU^,) + j dv2pU^i,vS)+ (40) 

^ j d^ndhi [Qa/3(ri, r2)/^^(ri2) - Pcpiji, ra) [fapijr 12)^12! ap^r 12) - U'ap{n2)) 



a,/3=p,n ' 

The solution of S < H2 > /SfafB = provides the Euler Correlations 
The quantities Q{P)ai3 are generalizations to the isospin dependent case of those of 
Ref.§: 

Qa,{r„ r2) = S^,^pUru r^) + ^Po (ri)p^(r2)V:,f (n^) - (41) 
Po,p(ri, r2)pj,p(ri, r2)V;°f^p(ri2) + Po,A(ri, r2)pj^(ri, r2)V;°f^^(ri2)] , 

Pa^i^u r2) = ^ {p^(ri)p^(r2) - 26^^ [(p^,p(r„ r2))^ + {plA^i, ^2)?] } , (42) 
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where: 



V:lp{r,2) = -25^f^[v^'\n2) + v^'Hn^) + Sv^'\n2) + ^v^'\ru)] - (44) 
4(l-Mb^'nn2) + 3t;(^)(ri2)], 

VZUin^) = -26^(^[v^'\n2) + v^'\n2) - v^'Hn^) - v^'\n2)] - (45) 

The contributions to < H2 > decouple in the (a/3) channels, so the minimizations 
are independently performed for each component. The corresponding Euler-Lagrange 
equations are: 

Kpi^) - Vafs{r)uafs{r) = Xaf3Uafs{r), (46) 

with: 



Ual3{r) 



1 



4P„;3(r) 



(47) 
(48) 



The Pa/3 and terms are obtained by integrating Pa^(ri,r2) and Qa/3(i'i5 1'2) over 
ri and r2 and keeping ri2 fixed: 



X(ri2) 



1 

ri2 Jo 



driTi 



|t'12-»'i| 



rfr2r2X(ri,r2,ri2) 



(49) 



where X can be either P or Q. 

Eqs. (|46|) are solved under the constraints: 

fap{r > dap) = 1 



'a/3 J 



0, 



(50) 



where the healing distances da/B are variational parameters. The same value for all the 
heahng distances {da/s = d) has been used in this paper. 

In some calculations, we have also adopted the Average Correlation Approximation 
(ACA), consisting in using an unique correlation, independent on the isospin of the nu- 
cleons, for both types of correlations, either gaussian or Euler. 

For the gaussian correlation, ACA simply consists in using the same values of the A 
and B constants for each channel. 

The ACA Euler equation is immediately obtained by summing over the a, (3 effective 
potentials. The resulting equation is: 



u"{r) — V{r)u{r) = Xu{r), 



(51) 



10 



where: 



u{r) = r^P{r)f{r), (52) 
1 



V{r 



AP{r) 



V'P{r) + 2Q{r) + ^ 
Fir] 



(53) 



Pir) = Y.Po.p{r), (54) 

a/3 

Q{r) = Y.Qap{r). (55) 

a/3 

The second ingredient of the calculations is the set of s.p. wave functions. In this 
paper, they have always been generated by a mean field potential of Woods-Saxon type: 



where m^r is the pion mass. 

As a test of accuracy of the calculations, we show in Tab. 2 the level of exhaustion of 
the previously discussed densities and distribution functions SR's. They have been calcu- 
lated in the FHNC/0 approximation, consisting in neglecting the elementary (or bridge) 
diagrams. These diagrams, designed as E^^ in Appendix A, cannot be summed in a closed 
way by FHNC Therefore, they are usually neglected. In some cases, their contribution is 
approximated by considering only a few low order terms. The seemingly crude FHNC/0 
approximation has been shown to be accurate in relatively low density systems as nuclei 
and nuclear matter P], whereas the elementary diagrams become increasingly important 
in high density systems as liquid atomic Helium [|TU]. On the other side, it has already 



been pointed out that the SR's are evaluated just to ascertain the degree of accuracy of 
the approximations employed in solving the FHNC equations. In the following, we shall 
always use FHNC/0, unless differently stated. 

The results shown in Tab. 2 have been obtained with the parameters of the ACA 
calculations of Tab. 5, whose details will be discussed later. In the non ACA case, the 
accuracy is substantially the same. The SR's are very well satisfied for all the medium- 
heavy nuclei. The worse situation is met in ^''^Pb, but the error remains less than 2%. 

In the last two rows of the Table we show the ratios between the values of the spin SR 
of Eq. (P^D obtained in the correlated system and that obtained in the uncorrelated one 
{S^""' / S'^'^^) . In the row labelled S'0-,0 the correlated value has been evaluated in FHNC/0, 
while in the row labelled S'0-,1 the contribution of the first order exchange elementary 
diagram has been included (see Ref.0 for a more extensive discussion of this point). This 
last approximation has been termed as FHNC-1. In agreement with the findings of Ref. 
0, Tab. 2 shows that the FHNC-1 approximation largely improves the accuracy of the 
calculated spin sum rules with respect to FHNC/0. 

A first set of calculations has been performed to investigate the relevance of the j- 
j coupling, of the Coulomb interaction and of the separated treatment of protons and 
neutrons. In order to have a better control on their contributions, we have used the 
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same set of Woods-Saxon parameters, identical for protons and neutrons and without 
spin-orbit term. In addition, we have taken the same s.p. potentials in ^^C and ^^O and 
for the Ca isotopes, as shown in Tab. 3, to investigate the influence of the unsaturated 1 
shells. The semi-realistic S3 interaction of Afnan and Tang [|l^ , supplemented in the odd 
channels with the repulsive interaction given by the repulsive terms of the even channels 
as discussed in Ref. [0, has been adopted in conjunction with the Euler ACA correlation 
function. 

In Tab. 4 we present the results of the ground state energies. Column Fl has been 
obtained by switching off the Coulomb interaction and the statistical correlation arising 
from the j-j coupling, Eq.(|l^). Its effect can be seen by comparing the results of column 
Fl with those of column F2, where the j-j coupling has been reinstated. Only ^^C, ^^Ca 
and 208 p]^ nuclei, which have some unsaturated 1 shell, are affected by this correlation and 
its influence turns out to be rather small. 

The F3 column shows the results obtained including the Coulomb interaction in the 
two-body hamiltonian. The contribution of the nuclear interaction {V in the Table) to 
the binding energy is about the same for '^°Ca, ^^Ca and ^''^Pb. These nuclei are large 
enough to allow the nuclear interaction to saturate. The contribution of the Coulomb 
interaction V^, because of the infinite range of the force, increases like the number of 
proton pairs, as expected. 

The results of the F4 column have been obtained by inserting the Coulomb potential 
also in the mean field. We have used the potential generated by an uniform spherical 
charge distribution. Other choices did not make any difference from the numerical point 
of view. 

A fully implemented variational principle would require performing the minimization 
of the energy functional respect to both the s.p. wave functions and the correlation 
function parameters. This implies a large parameter space minimization. Because the 
requested numerical effort is heavy, and also because we are not yet considering fully 
realistic interactions, the set of single particle wave functions has been kept fixed and the 
minimization has been done only respect to the correlation function. 

Moreover, for all the columns of Tab. 4 we have used the same correlations (obtained 
by minimizing the Fl energies). Therefore, we cannot draw any definite conclusion on 
whether the 1-s coupling is or is not variationally preferred with respect to the j-j coupling, 
except that we find very small differences in the two cases. Similarly, we should not be 
surprised by the fact that the inclusion of the Coulomb effects in the wave function (column 
F4) does not seem to lower the energies with respect to their perturbative estimates of 
column F3. 

In all the calculations we shall discuss henceforth the s.p. wave functions have been 
generated by Wood-Saxon potentials whose parameters, given in Tab. 5, are taken from 
literature 0, and have been fixed to reproduce the s.p. energies around the Fermi surface 
and the rms charge radii. 

Table 6 gives the results obtained with the S3 interaction and various correlation 
functions. The EU rows show the results with the isospin dependent Euler correlation 
function, the ACA rows those with the ACA Euler correlation and the G rows those with 
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the gaussian ACA correlation. The columns present the contributions of the various terms 
building up the average value of the binding energy. The binding energies per nucleon, 
E/A, have been calculated by subtracting the center of mass energy. In the last column 
the healing distances, for the minimum of the energy with the Euler correlations, are 
given. Note that the EU calculations have been performed using the same value of the 
healing distance for all the isospin channels. 

For all the nuclei we have considered, the isospin dependent correlation function re- 
sults to be variationally preferred respect to the simpler ACA functions. This is mainly 
produced by a large decrease of the potential energy, which overcompensates the increase 
of the kinetic energy. The poorest choice is the ACA gaussian one. 

The correlation functions used in the ACA Euler and gaussian calculations of Tab. 6 
are shown in Fig.l. They are very similar in all the nuclei we have studied and all the 
finite system Euler correlations overshoot one. 

Also in the isospin dependent case, where /^^, /"'^ and /"^ are different even if we 
assume the same healing distance for all of them, the correlations do not show a strong 
dependence on the type of nucleus considered. This is due to the fact that the correla- 
tions are mostly sensitive to the short range behaviour of the nuclear interaction. Fig. 2 
compares the isospin dependent Euler correlations with the ACA Euler ones for ^^C and 
^°^Pb. /^^ and are close to f^'~^^i but f^'^ shows large differences. This rich isospin 
structure is responsible for the better variational behaviour of the correlation. 

It is interesting to compare the finite nuclei Euler correlations with those of infinite 
nuclear matter. In this case, we put = /^^ = f'^j} and Z"'' = \ [3/^^/ + fJjM] 5 where 
T = 0, 1 is the total isospin. The nuclear matter correlations are always smaller than 
one, therefore the overshooting feature of the Euler correlations in the nuclei seems to 
be ascribed to the finite nature of the systems. Nuclear matter correlations have been 
often used in finite nuclei [0 on the assumption that their short range behaviour is 
little affected by the size of the system. However, we find not negligible differences and 
improved energies by using correlations generated by the finite system Euler equations. 
For instance, in ^^C /^a/ gives E/A = —2.12MeV at d = 2.0/m, to be compared with 
the E/A = —4:.65MeV value of the Euler correlation. 

The proton one-body densities, Pi(r), calculated with the Euler and gaussian cor- 
relations are given in Fig. 3 for ^^C, ^^O and the Calcium isotopes and Fig. 5 for ^°^Pb, 
and compared with their uncorrelated counterparts, Po(r), (full lines). The dotted curves 
have been obtained with the gaussian-ACA correlations and all of them are lower than 
the uncorrelated densities at low r-values. On the other hand, those obtained in the Euler 
(dashed-dotted lines) and ACA Euler (dashed lines) calculations are larger than Pq^t). 
This indicates that the correlation effect on Pi(r) at small distances may strongly depend 
on the type of correlation employed. The results of the figures have been obtained for 
those values of the variational parameters giving the minimum energy. We found corre- 
lated densities smaller than Po{r) in the nuclear center when we used values of the healing 
distance far from those producing the energy minima. 

Similar results have been obtained for the neutrons density distributions. 

In order to check the validity of our approximations and their sensitivity to the choice 
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of the nucleon-nucleon interaction, we have performed ACA Euler calculations with the 
Brink and Boeker Bl central potential fl^. The Bl potential is not a microscopic one, as 



it does not fit two-body data. It is an effective interaction which reproduces the nuclear 
matter and ^He binding energies in Hartree-Fock theory. This potential has been used 
in Ref. ||^ to compare its FHNC results with the available, exact Variational Monte Carlo 
calculations in ^^O. 

The accuracy of the SR's for the Bl potential is the same as for S3. The ground state 
expectation values of the various terms of the hamiltonian are shown in Tab. 7. Also 
in Bl case, the correlation functions, shown in Fig. 6, are very similar for all the nuclei 
considered. However, some differences from those obtained with the S3 interaction are 
present. The Bl correlations have a value of about 0.55 at r = fm, versus the 0.37 value 
of the S3 case. In addition, the overshooting is rather small. 

The differences in the correlation functions have some consequences on the proton 
density distributions, shown in Figs. 3 and 5 by the dashed-doubly-dotted lines. All these 
distributions are smaller than the uncorrelated ones, in contrast with the ACA Euler 
results with S3. 

We have studied the effect on the energy of the insertion of the dressed, lowest order 
exchange elementary diagrams in FHNC-1. The results are summarized in Tab. 8, where 
the FHNC-1 energies per nucleon (columns Ei) obtained with the Euler-ACA correlations 
for both S3 and Bl interactions are given. In addition we give the percentile deviations 
A from the FHNC/0 energies: 

A = 100^^5^ (57) 

where Ei and Eq are the FHNC-1 and FHNC/0 energies respectively. 

Confirming the finding of Ref.|p, Ei is always higher than Eq. Moreover, the influence 
of this class of elementary diagram becomes small in large nuclei, in agreement with 
nuclear matter results |T^. 

The proton MDs, n^{k), obtained for S3 are shown in Fig.4 for ^^C,^^0,^°Ca and 
^^Ca nuclei, and, in the lower panel of Fig. 5 for ^''^Pb. In the same figure we show the 
uncorrelated distributions (full lines) and nPlk) for the Bl model with the ACA Euler 
correlations. The MD sum rules of Eqs. (|37|j38| ) are given in Tab. (9) for different wave 
functions. The MDq sum rules are very well satisfied, with errors less than 1%, except 
one single case. The MD2 sum rules are also well satisfied, even if with a lower accuracy. 
However, for the best variational case (the Euler correlation) the error in MD2 is, at most, 
about 5%. 

All the FHNC results show the same known features: the correlated results are orders 
of magnitude larger that the uncorrelated ones at high momentum values. In particular, 
the isospin dependent correlations enhance the tails of the MDs more than the Jastrow, 
isospin independent choice. This result is in line with the findings of Ref. 0, where the 
effect of state dependent correlations (having spin, isospin and tensor components) was 
estimated, from CBF Nuclear Matter calculations, in a suitable Local Density Approx- 
imation and compared with a Jastrow finite nuclei FHNC evaluation of n{k) in N=Z 
nuclei. State dependence was found to increase the tails of the MDs of about one order 
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of magnitude more than the Jastrow correlations. The enhancement found here is some- 
what smaller. However, a consistent FHNC approach, to deal with nuclear matter type 
correlations, is called for before a quantitative an definite answer may be given. 

5 Conclusions 

This work represents a further step towards a fully microscopical description of medium- 
heavy nuclei. For the first time we have performed FHNC calculations for medium and 
heavy finite nuclear systems. The calculations have been done using different s.p. wave 
functions for protons and neutrons, in j-j coupling, isospin dependent correlations and 
semirealistic, central potentials. The Coulomb interaction has been also considered. 

The extension of the FHNC equations of Refs.]^, P] is not straightforward because, 
due to the j-j coupling scheme, a new type of statistical link appears. 

We have used the central S3, Afnan and Tang, and Bl, Brink and Boeker, interac- 
tions. The s.p. wave functions have been generated by Woods-Saxon potentials, whose 
parameters have been taken from literature. 

A detailed analysis of the influence of the j-j statistical correlations on the FHNC 
equations has been presented. Their contribution is rather small, while the effect of the 
Coulomb interaction, especially for heavy nuclei, turns out to be relevant. 

Both isospin dependent and independent correlation functions have been used. The 
first ones allow us to distinguish between pp nn and np channels. In all the nuclei consid- 
ered we have found that the pp and nn correlations are quite close to each other, while 
the np correlation shows pronounced differences. This marked isospin structure appears 
to be variationally preferred, as it largely lowers the ground state energies with respect 
to the results obtained with isospin independent correlations. 

At the variational minima, the correlation functions do not strongly depend on the 
considered nucleus. However, they show an intermediate range structure, not present in 
the optimal FHNC correlations of infinite, symmetric nuclear matter. This difference has 
large effects on the energies and the density distributions. 

The accuracy of the FHNC/0 and FHNC-1 approximations has been studied by exam- 
ining the density and distribution functions sum rules. FHNC/0 is found rather accurate, 
the largest error being less than 2% in the scalar sum rules. The inclusion of the lowest 
order, exchange elementary diagram, in FHNC-1, lowers the error in the spin sum rules 
to less than 1%. The FHNC-1 correction is also relevant in the ground state energy, 
because the adopted potentials have large Majorana parts [Q. However, its contribution 
is negligible for large mass nuclei, as ^°^Pb. 

As far as the matter density is concerned, it appears that, given a chosen s.p. structure, 
the correlations may largely and differently modify the density distribution. The results we 
have presented are however related to energy minima obtained by performing variations of 
the correlation functions only. Before drawing any sensible conclusion one should include 
tha correlations in a fully consistent scheme. In our case, thise means to perform the 
energy minimization by variating also the mean field parameters. 
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The momentum distributions have been calculated within the FHNC/0 scheme and 
the sum rules analysis shows that its accuracy is actually very satisfactory. Short range 
correlations produce high momentum tails, absent in an independent particle model. The 
momentum distributions calculated with state dependent correlations show tails higher 
than those produced by Jastrow, state dependent correlations. 

In this paper, we have shown, that it is possible to treat accurately doubly closed shell 
nuclei in the j-j scheme using the FHNC cluster summation technique. This makes feasible 
the microscopical investigation of heavy nuclei, such as ^'^^Pb. Isospin state dependence 
in the correlation can be safely handled and provides a better variational choice than 
Jastrow correlated wave functions. This encourages us to pursue the goal of a microscopic 
description of finite nuclei with even more realistic correlation functions and interactions 
containing, in particular, tensor components. 
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Appendix A 

This appendix contains the two-body distribution function and one-body density FHNC 
equations for the j-j couphng and isospin dependence. The topological classification of 
the cluster diagrams (Nodal, Composite and Elementary, dd, de, ed and cc) has been 
already presented in Ref. 0, as well as the chaining technique; so it will not be repeated 
here. 

In the following we shall indicate the coordinate with i, and the convolution integral 
over the coordinate 3 as (|). 

The sums of the nodal diagrams N'^y=dd,de,ed,ee ^^re given by the solutions of the integral 
equations: 

A^.f (1, 2) = E E 3)C\'(3) I iV/,(3, 2) + X/^(3, 2)) , 

7=P>"a;',j/'=d,e 

where = otherwise Qy = Q. The sum over {x'y') runs over the pairs {dd,de,ed), 
because of the FHNC convolution rules. 

The cc-chains now can contain the parallel exchange link p{l, 2) and the antiparallel 
one Pq 2). We give here their explicit expressions, as they represent the main changes, 
due to the j-j coupling, respect to the usual symmetric matter FHNC Eqs. 

iVi:);(l, 2) = 3)C(3) I AT,", p(3, 2) + X:,,p(3, 2) - p^,p(3, 2)) - 

{X:^JI, 3)C(3) I 2) + X:^J3, 2) - p°^(3, 2)) 

iVi$(l, 2) = - (p?,p(l, 3)C(3) I <^;(3, 2) + X,",,p(3, 2)) 

- (Po,p(l, 3)(C(3) - 1) I iV^,^;(3, 2) - p^,p(3, 2)) 
+ {pl^l, 3)C(3) I <!4"(3, 2) + X:,J3, 2)) 

+ (pm(1'3)(C(3)-1) I <);(3,2)-p°^(3,2)) 
2) = {X:,J1, 3)C(3) I NZ^piS, 2) + X:,,p(3, 2) - p^,p(3, 2)) + 
(X,",,p(l, 3)C(3) I N:,J3, 2) + X,V(3, 2) - p^,^(3, 2)) 
ivS(l, 2) = - (p?,^(l, 3)C(3) I (3, 2) + X,",,p(3, 2)) 

- (po,a(1, 3)(C(3) - 1) I iVi^J (3, 2) - p^,p(3, 2)) 

- (Po,p(l, 3)C(3) I <lr(3, 2) + X,",,^(3, 2)) 

- {PoAh 3)(C(3) - 1) I <5(3, 2) - p^,^(3, 2)) , 

and 

iV"c,x(l, 2) = iVi:S(l, 2) + iVi^5(l, 2). 
The partial pair FHNC distribution functions, g^y, and the sums of the composite 
diagrams, X^y, are: 

^?f(l,2) = /i(l,2)exp 



2) +i?2f (1,2) 
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= 1 + <(1,2)+X,7(1,2) 



a/3/ 



5f(l,2) = ^,7(2,1) 



rf(l,2)[Arf(l,2) + £;,t(l,2) 
iV.f(l,2) + Xf(l,2) 



^^^(1,2) = ^f(l,2)[iV,f(l,2) + ^f(l,2) + 

(TV^f (1, 2) + E:S{1. 2)) 2) + (1, 2)) 

25„;3 (Ar.V(l, 2) + 2) - P^,p(l, 2))' - 

2<^a^ (a^cV(1' 2) + ^cV(l' 2) - PoV(l, 2) 
7V«/(1,2) + X,f(l,2) 



Cx(l,2) = ^,7(l,2)[A^eV(l'2) + £;^,x(l,2)-PoV(l'2) 
= 7V-^(l,2)+X,",,^(l,2)-p^,;,(l,2). 

The vertex corrections, Ci=d,e> given by: 

= exp 

= C(l)[t^"(l)+Po(l)], 

where 



/ ^'^2 E {e.'(2) [^2f (1, 2) - (1, 2) - 5f (1, 2) (^f (1, 2) - l)" 

2) - 2) - 5f (1, 2) 2) - l) 



+ 



ef(2) 



^dt(l,2)^r(l,2) 



a/3/ 



+ E2{1) 



Ut{l) = / ciVJ E |e.'(2) kf (1, 2) - E:^{1, 2) - 2) 2) - l) + 



+ 



ef(2) 



(1, 2) - (1, 2) - (1, 2) 2) - l) 

^et (1, 2)^?f (1, 2) - (1, 2)5ef (1, 2) - 



'5dt(l,2Kf(l,2) 



a/3/ 



+ 
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4C(2) 



^1p(1, 2)Cp(1, 2) + S:,J1, 2)g:Jl, 2) 



P^,p(l,2)(<V"(l,2)-p^,p(l,2)) + 



Po,a(1,2)(<);(1,2)-p"^(1,2) 



with: 



2) = 2) + (l/2)iV"/^(l, 2 



Notice that p°(l) = ^^(1) 



Appendix B 



Here we will shortly present the FHNC equations for the one-body density matrices 
p''*"(l,l'). These are modifications of the equations presented in Ref.0. As we have 
discussed throught the text we have to separate the case with the same third component 
of spin 

P^^"(i,i') = -C(i)C(i')^?l(i,i')[A^"....p(i,iO +K...,p(i,i') -Po.p(i,i')], 

from the case with different ones 



p--(l,l') 



g"^ is defined as: 



exp 



The nodal functions N^^, N^^^^ x solutions of convolution equations, formally iden- 
tical to those presented in appendix A, with the substitution of the external (1, 2) link 
indices {d,c) by {lj,uJc), where needed. New nodal terms N"d uje^ ^!^cc x appear in these 
convolutions. The same rule as before applies to them. 

New partial FHNC distribution functions need to be introduced. They are given by: 



£."(1,2) 



^?i"(2,l) 
/„/3(l,2)exp 

1 + Cf(l,2) 



Cf(l,2)+Kf(l,2) 



^.1(1,2) 



£f(l,2) 



/3q 
af3 



^7ej(2,l) 



£."(1,2) 



Ar"fri,2) + R"fri,2i 



Ar^f(l,2) + Xl^(l,2) 



£rc,x(l,2) 



£."(1,2) 



AT- ^(1, 2) + EZ,xil, 2) - p^,^(l, 2) 



AT^"" ^(1, 2) + XZA^^ 2) - p" (1, 2) 
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The vertex corrections, is: 



C(l) = exp[t/:(1)], 



with: 



0=p,n 



Ki^) = / dh, Y: [^,{2) [xfA^, 2) - EfA^, 2) - SfAh 2) (^^.^(l, 2) - l)" 

2) - (1, 2) - (1, 2) (£,^(1, 2) - l) 



+ 



ef(2) 



Appendix C 



In this appendix, we shall give the explicit expressions for the uncorrelated one-body 
densities and related quantities. For the s.p. wavefunctions we use: 

IJ,,s ^ 



The one-body density pj.^ is given by: 
1 



P?i(ri) = ^E(2j + l) 

nlj 

where we have defined: 



/\2 



'1 '1 
For the one-body density matrix we find: 

PT2(ri,r2) 



2(47r)2 



^(2j + l)(2/ + l)i?;j,^.(r2)i?:,,,.,(r2) 



n'l'j' 



^<;,(riXT,^(ri)P/(cos^)i^Hcos^)| + 



(47r) „y 

n'l'j' 

I [<,,.(ri)D;j,,^.,(ri) - i?;^;,.(riX,V,.,(ri)] gKcos^)gr(cos^) 
sin2(^) 



-Kij{n)Kn'A^i)Q'i{cos9)Q'Acos9)\, 
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where Pi are Legendre polynomials, 9 is the angle between the vectors ri and r2 and we 
have defined: 

Q/(cos^) = sin^P/(cos^) 

q;(cos^) = ^(cos^P/(cos^) -/(/ + l)Pz(cos^)). 
sm u 

For the densities we obtain: 



nlj 



nlj 



Finally, pj.4, appearing in the calculation of the center of mass term Tcm, is given by: 
PT4(ri, r2) = PT6(ri, ^^2) - Po,p(ri, V2)Pt^^p{yi, V2) - Po,^(ri, r2)pr5_^(ri, ra), 

where 

PT6(ri, r2) = 2 (ViPo,p(ri, • V2Po,p(ri, + ViPo,^(ri, r2) • V2Po,A(ri, r2)) , 
PT5,x(ri,r2) = 2Vi • V2Po,x(ri,r2), 



and 



nlj 



^RUr,)?^ + R^n)?^] sin^^P/(cos^) + 



ri 



r2 



Rnij{ri)Kij{r2) 



27r 



rir2 

E(-i) 

nlj 

nlj ( 



(sin2^P/(cos^) + /(/ + 1) cos^Pi(cos^ 



j-i-l/2 



^:0(n)^:z',(^2)cOS^gKcOS^) + 



<z,(r2)^?^ + sin^^QKcos^) 



ri 



r2 



<^ sm 9Q, (cos 6')+ 

rir2 >- " 

)S 



Z(Z + 1) 



sin"^ 



PT6(ri,r2) = 



2(47r)2 



E(2j + l)(2/ + l)P;:,,.(r2)P:,,,.,(rO 
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cos^ 



Kij ir2)Ki'j' iri)Pi (cos 9)Pi, (cos 6) - 

+ 



• 2 Q^nlj{f'i)RZ'l'3'{fl) -r,l( Q\ dU 

sm Pi (cos6')P;,(cos 



sin^^ 



Rnljir2) 



rir2 

?^^R::i,^,{ri)Pl{cos9)Pu(cose) + 



Pi{cos9)P{,{cose) 



(47r)2 % 



T.i-^y^''-'-'-'Rnijir2)Ki'A^i) 



COS 9 



KUr2)R::,An)fi{cos9)Mcos9)- 



■ 2a^rilj{^2)Rn'l'j'{^l)fi, a\ fi t o\ 

sm^ 9 — fi (cos 9)fi,{cos 9) 



rir2 



+ 



sin'^ 



Rnljir2) 



Rnlj{r2) 



r2 



R^n'A^i)n{cos9)Mcos9) + 



fi{cos9)f[,{cos9) 
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Table 1. Coefficients of eq. (ID . 





12Q 


WQ 


40Ca 


48Ca 


208pj^ 


Sp 


0.999 


1.000 


1.000 


1.000 


0.995 




0.999 


1.000 


1.000 


0.999 


0.996 


>Jpp 


0.994 


0.998 


1.000 


0.999 


0.989 


>-'np 


0.998 


1.000 


1.000 


0.999 


0.991 


Snn 


0.994 


0.997 


1.000 


0.998 


0.993 




1.004 


1.018 


1.009 


1.014 


0.953 


Scr,l 


1.001 


1.000 


1.000 


1.002 


0.985 



Table 2. Ground state sum rules. ^^,o{i) are the ratios S^f^'^'/S^'"' in FHNC/0 and 
FHNC-1 approximations, respectively (see text). 
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R 


a 




52.5 


3.20 


0.53 




52.5 


3.20 


0.53 




57.5 


4.10 


0.53 




57.5 


4.10 


0.53 


208 


60.4 


7.46 


0.79 



Table 3. Coefficients of the Woods-Saxon potential (|56D used for the calculations of 
Tab.4. 
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12^ 


Fl 


F2 


F3 


F4 


V 


-21.28 


-21.43 


-21.43 


-21.24 


Vc 


0.00 


0.00 


0.64 


0.63 


T 


19.10 


19.28 


19.28 


19.09 


E 


-2.18 


-2.15 


-1.51 


-1.52 


160 


Fl 


F2 


F3 


F4 


V 


-26.18 


-26.18 


-26.18 


-25.84 


Vc 


0.00 


0.00 


0.87 


0.86 


T 


20.16 


20.16 


20.16 


19.87 


E 


-6.02 


-6.02 


-5.15 


-5.11 




Fl 


F2 


F3 


F4 


V 


-32.61 


-32.61 


-32.61 


-31.97 


Vc 


0.00 


0.00 


1.95 


1.91 


T 


23.96 


23.96 


23.96 


23.45 


E 


-8.65 


-8.65 


-6.70 


-6.61 




Fl 


F2 


F3 


F4 


V 


-33.41 


-33.69 


-33.69 


-32.13 


Vc 


0.00 


0.00 


1.62 


1.59 


T 


25.77 


26.12 


26.12 


25.68 


E 


-7.64 


-7.57 


-5.95 


-5.86 


208 


Fl 


F2 


F3 


F4 


V 


-33.44 


-33.70 


-33.70 


-32.67 


Vc 


0.00 


0.00 


3.98 


3.84 


T 


24.04 


24.29 


24.29 


23.68 


E 


-9.40 


-9.41 


-5.43 


-5.15 



Table 4. Energies per nucleon, in MeV , for the five nuclei considered. These results 
have been obtained using the Woods-Saxon parameters of Tab. 3, the S3 interaction of 
Afnan and Tang []rT| and the Euler ACA correlation. The healing distances are d = 2.1 fm 
for ^^C and d = 2.0/m for the other nuclei. 
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H Vls R a 




P 
n 


62.00 3.20 2.86 0.57 
60.00 3.15 2.86 0.57 


160 


P 
n 


52.50 7.00 3.20 0.53 
52.50 6.54 3.20 0.53 




P 
n 


57.50 11.11 4.10 0.53 
55.00 8.50 4.10 0.53 




P 
n 


59.50 8.55 4.36 0.53 
50.00 7.74 4.36 0.53 


208pb 


P 
n 


60.40 6.75 7.46 0.79 
44.32 6.08 7.46 0.66 



Table 5. Coefficients of the Woods-Saxon potential ([56|). 
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-335.5 


8.2 


228.7 


65.4 


12.8 


-33.2 


-3.84 


2.44 
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-291.9 
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45.0 
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-15.5 


-2.36 
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273.3 


82.1 


11.4 


-119.7 


-8.20 


2.36 


160 


ACA 


-418.5 


13.9 


266.8 


55.0 


11.4 


-82.8 


-5.89 


2.08 




G 


-395.8 


13.5 


263.2 


43.4 


11.4 


-75.7 


-5.44 






EU 


-1515.1 


80.4 


791.3 


262.6 


10.1 


-381.2 


-9.78 


2.42 


40Ca 


AGA 


-1302.2 


78.0 


771.9 


190.2 


10.1 


-262.1 


-6.81 


2.20 




G 


-1220.0 


75.6 


750.3 


154.2 


10.1 


-239.9 


-6.25 






EU 


-1744.9 


78.8 


968.8 


302.4 


9.6 


-394.9 


-8.43 


2.45 


^^Ca 


AGA 


-1494.5 


76.0 


936.6 


209.6 


9.6 


-272.3 


-5.87 


2.15 




G 


-1418.8 


74.3 


918.7 


182.7 


9.6 


-243.2 


-5.27 






EU 


-7600.3 


810.8 


3865.3 


1163.6 


6.0 


-1761.4 


-8.50 


2.52 


208pb 


AGA 


-6772.0 


809.2 


3910.7 


995.8 


6.0 


-1056.3 


-5.11 


2.45 




G 


-6218.4 


792.1 


3740.7 


814.1 


6.0 


-871.5 


-4.22 





Table 6. Ground state expectation values with different correlations. EU indicates the 
Euler correlation, ACA the ACA Euler correlation and G the ACA gaussian correlation. 
The parameters of the gaussian correlations are A = 0.7 for all the nuclei considered and 
B = 2.15/m-2 for ^^Ca, B = 2.2fm-^ for ^^C, ^OQa and ^ospb, and B = 2.3/m-2 for 
-"^^O. = t9'^ + tP-* + Tp-*. Energies in MeV and distances in fm. 
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< y > 


<v,> 




Tf 


T 

cm 


< H > 


E/A 


d 


1 o ^ < 
12Q 


-293.6 


8.1 


221.9 


15.6 


12.8 


-48.0 


-5.07 


2.17 


160 


-435.0 


13.6 


252.4 


20.9 


11.4 


-136.7 


-9.26 


2.07 




-1328.6 


76.7 


759.2 


78.2 


10.1 


-414.5 


-10.62 


2.25 




-1554.1 


75.6 


931.1 


93.2 


9.6 


-454.2 


-9.66 


2.30 




-7017.8 


805.3 


3827.5 


390.2 


6.0 


-1994.8 


-9.62 


2.50 



Table 7. Ground state energies for the ACA Euler correlation and the Bl potential. 
Energies in MeV and distances in fm. 





S3 


Bl 




El A 


El A 


120 
160 

40Ca 

48 Ca 

208p|j 


-2.21 6.36 
-5.57 5.43 
-6.33 7.04 

-5.28 10.05 
-5.14 -0.59 


-4.63 8.68 
-8.61 7.02 
-10.05 5.37 

-9.33 3.42 
-9.69 -0.73 



Table 8. FHNC-1 binding energies per nucleon calculated with the Euler-ACA corre- 
lations for the S3 and Bl interactions. A is the percentile deviation from the FHNC/0 
energy. 
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MD2/ <T> 




EU 


0.990 


0.990 


0.955 




ACA 


0.995 


0.995 


0.960 




G 


0.997 


0.997 


0.942 




EU 


0.996 


0.996 


0.969 


160 


ACA 


0.996 


0.996 


0.963 




G 


0.996 


0.996 


0.954 




EU 


1.000 


0.995 


0.956 




ACA 


0.995 


0.996 


0.942 




G 


0.992 


0.993 


0.919 




EU 


0.998 


0.993 


0.949 


48Ca 


ACA 


0.993 


0.994 


0.936 




G 


0.991 


0.993 


0.914 




EU 


0.999 


1.005 


0.943 




ACA 


0.996 


1.002 


0.934 




G 


0.984 


0.985 


0.917 



38|), for the calculation of Tab. 6. Xp 



Table 9. Momentum distributions SRs, Eqs. 
Z/A and x„ = N/A are the proton and neutron concentrations respectively. 
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Figure captions 



Fig.l. ACA Euler (upper panel) and gaussian (lower panel) correlation functions at the 
energy minimum for several nuclei, for the S3 potential. 

Fig. 2. Euler and ACA Euler correlation functions at the energy minimum for ^^C and 
^®Ca nuclei with the S3 interaction. Full lines are the ACA Euler correlations. The lowest 
panel gives the nuclear matter nn{— pp) (dot-dashed) and np (dashed) Euler correlations 
fot the same potential, at d = 2.0/m. 

Fig.3. Proton density distributions. Full lines are the uncorrelated densities. The dotted 
curves are obtained with the ACA gaussian correlation, the dashed curves with the ACA 
Euler correlation and the dot-dashed curves with the isospin dependent Euler ones. The 
S3 interaction has been used. The dash-doubly-dotted curves are the densities obtained 
with the Bl potential using the ACA Euler correlation. 

Fig.4. Proton momentum distributions. Curves as in Fig.3. 

Fig. 5. Proton density (upper panel) and momentum distribution (lower panel) for 
^°^Pb. Curves as in Fig.3. 

Fig. 6. ACA Euler correlation functions for the Bl interaction. 
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r [fin] 




r[fm] 




0.0 1.0 2.0 3.0 4.0 




0.0 1.0 2.0 3.0 4.0 

k [fm'] 




